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Abstract 

Stochastic and bistochastic matrices providing positive maps for spin 
^ ' states (for qudits) are shown to form semigroups with dense intersection 

^\ ' with the Lie groups IGL{n, R) and GL{n, R) respectively. The density 

CNj , matrix of a qudit state is shown to be described by a spin tomogram de- 

CO ' termined by an orbit of the bistochastic semigroup acting on a simplex. A 

^^ I class of positive maps acting transitively on quantum states is introduced 

r~~^ . by relating stochastic and quantum stochastic maps in the tomographic 

("^ ' setting. Finally, the entangled states of two qubits and Bell inequalities 

OO , are given in the framework of the tomographic probability representation 

^^ ■ using the stochastic semigroup properties. 

Key words Quantum tomograms, semigroups, positive maps, entangled 

states. 
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1 Introduction 

The description of a physical system admitting a probabilistic interpretation, 
be it classical or quantum, requires two collections of objects called states and 
observables, say S and O respectively, along with a pairing jj. associating with 
any state p and observable A a Borel probability measure on the real line M. If 
A is measured while the system is in a state p, fJ.A,p represents the probability 
distribution for the observed values of A. Thus if i? C R is a Borel set, fiA.p{E) S 
[0, 1] is the probability that the measured value of A will be in the set E when 
the system is known to be in the state p. From a general point of view, above 



properties seem to be the minimal features that any physical system should 
possess. 

This approach has been studied by several authors, for instance one can find 
a nice discussion by G. Mackey[T]. The set S describes the basic mathematical 
structure we are dealing with, while fiA,p provides us with a physical interpre- 
tation. This probabilistic point of view is compatible with convex combinations 
on the space of states, indeed if pi and p2 give rise to probability distributions, 
by setting ^J-A,\pl+{l~^}p2 = ^MA.pi + (1 - '^)ma,p2 we define a new probability 
distribution when < A < 1. Usually one requires some additional structure 
telling us how the system changes from time s to a later time t, i.e. requires the 
existence of a family of mappings Ut^s '■ 5 — > 5 representing the dynamics and 
called evolution operator. The requirement that a state at a given time deter- 
mines the state at a later time forces us to postulate the semigroup property 

Ut,.t, = Ut,,s o U,.t, (1) 

with Ut^t the identity. Within this setting 

ji : S y. O ^ {Borel probability measures on M} (2) 

A subset of observables is said to be a tomographic set r if it allows to identify 
the state p (to "reconstruct" the state) when {/x^^pj^g^ is known. Very often r 
is generated by acting with a group Q on some fiducial observable Aq^ i.e. it is 
the orbit of t/ in O through Aq. According to the group we use and the fiducial 
observable we start with, we deal with symplectic tomography, photon-number 
tomography and so on. While this approach is general enough to allow us to 
deal both with classical and quantum systems, here, to be more definite, we 
shall consider a quantum system with a finite number of levels. 

States for quantum systems with a finite number of levels will be thought of 
as the spin states (or qudits) they can be described by density matrices which 
are hermitian nonnegative (2j + l) x (2j + l) matrices with unit trace. The linear 
maps of the spin states, positive maps, can be described by (2j -I- 1)^ x (2j -I- 1)^ 
matrices with special properties [5]. Recently it was shown [31 SI [S] that qudit 
states can be described by probability distributions of random spin projection 
(called tomogram) depending on the direction of the quantization axis. In view 
of this the geometry of qudit states can be associated with the geometry of 
a simplex and the set of positive maps of qudit states can be associated with 
stochastic and bistochastic matrices moving points on the simplex. The aim of 
this work is to find the connection of spin tomograms with a unitary matrix 
containing eigenvectors of the density matrix of a qudit state and a point on 
the simplex which has the eigenvalues of the density matrix as its coordinates. 

Another aim of the work is to define positive maps of qudit states through the 
transitive actions of both the unitary group on the eigenvectors of the density 
matrix and the stochastic matrix semigroup on the eigenvalues of the density 
matrix regarded as points of the simplex. 

The qudit states of multipartite systems can be either separable or entangled. 
We formulate the properties of a qudit tomogram, which is the joint probability 



distribution of two spin projections on their own quantization axes, able to dis- 
tinguish separable and entangled states. We consider the Bell inequalities [6l [7] 
in the context of the properties of stochastic matrices constructed by using spin 
tomograms. The Cirelson[5] bound 2\/2 for the Bell-CHSH inequahty of two 
qubits will be connected with some properties of a universal stochastic matrix 
obtained from the tomographic probability distribution describing maximally 
entangled two spin-1/2 states. The connection of positive maps with the semi- 
group of stochastic matrices provides the possibility to find a new relation of the 
maps with the Lie group of the general linear real transformations GL{n, M.) for 
bistochastic matrices and with the inhomogcncous group IGL{n,M.) for stochas- 
tic matrices. 

This connection (which seems to have been unknown) provides a possibility 
to construct unitary representations of stochastic and bistochastic semigroups 
by reducing known infinite dimensional unitary irreducible representations of 
the Lie groups to the subsets of the Lie groups which are the semigroups under 
consideration. 

The paper is organized as follows: in section 2 we review the spin tomography 
approach for one and two qudits. Examples of a qutrit and two qubit states in 
tomographic probability representation are studied in section 3. The relation of 
stochastic and bistochastic semigroups with Lie groups is discussed in section 
4, mainly in the case a qutrit. In section 5 a class of positive maps acting 
transitively on quantum states is introduced by relating stochastic and quantum 
stochastic maps in the tomographic setting. The relation of stochastic matrices 
with Bell inequality violation for entangled states of two qubits is discussed in 
section 6. Some conclusions and perspectives are finally drawn in section 7. 

2 Spin tomograms and unitary group 

As it was shown in [3l [9l [10] the qudit state described by a (2j -I- 1) x {2j + 
1)— matrix p can be also described by a tomographic probability distribution 
function, or tomogram, yv{m, U) > where m is the spin projection: m = 
— j, — j + 1, . . . , j — 1, j; and t/ is a unitary (2j -|- 1) x {2j + 1)— matrix. This 
matrix can be considered as a matrix of an irreducible representation of the 
rotation group depending on two Euler angles 0, 9 determining the direction of 
quantization (or a point on the Bloch sphere S'^). The physical meaning of the 
tomogram yV{m, U) is that, in the spin state with the given density matrix p, it 
gives the probability to obtain m as spin projection on the direction determined 
by the two angles 4>, 9. It corresponds to choose {[/'I' J^C/} as tomographic set of 

isospectral observables, where J^ = X]m=-i '^ 1"^) ("^1 i^ '^'^'^ of the generators of 
the irreducible representation of the rotation group, so that ^{7x1^ U) is nothing 
but the value of the concentrated measure p-jji j^u,p ^^ the spectral point m : 

Wim, U) = Hufj^uA^) = TrC/^ \m) {m\ Up = (m |C/pC/^| m) . (3) 



The probability distribution is obviously nonnegative and normalized, i.e. 

i 
^ >V(m, C/) = 1 (4) 

for any direction of the quantization axis. The spin tomogram can be also 
regarded as the diagonal matrix clement of the rotated density matrix UpU^ in 
the natural basis \m) . 

The relation is invertiblc and knowing the tomogram yV{m, U) for the ma- 
trices U{4), 6) of an irreducible representation of SU{2) one obtains the density 
matrix p by means of a linear transform [51 110| which is the analog of the integral 
Radon transform but in the space of qudit states. Thus the quantum state of 
a a qudit (a spin-j state) is known if the probability distribution W{m^ U) of 
random spin projection as a function of the unitary matrix U is known. The 
tomogram yi>{m, U) can be used, consequently, in alternative to spinors (wave 
functions) or density matrices for describing spin states. The information on 
the spin state contained in the tomogram is redundant since it is sufficient to 
know the tomogram only for several directions determined by a set of angles 
{4'k^0k\^ whose number corresponds to the number of parameters determining 
the density matrix, equal to (2j + 1)^ — 1. But at the same time the dependence 
of the tomogram yV(?n, U) on the parameters of the unitary matrix U provides 
some advantage in considering the spins, j = 0, 1/2, 1, . . . ; and also the quantum 
states of several spins in an unified approach. For two qudits (spin ji and J2) 
the tomogram of the quantum state with the (2ji + 1) (2J2 + 1) x (2ji + 1) (2J2 + 1) 
density matrix p is the normalized joint probability distribution 

W(TOi,m2, C/) = (mim2 |t/pC/^| mim2) (5) 

of two random spin projections mi ~ —ji, —ji + 1, ■ ■ ■ , ji ~ l,ii and 7712 = 
—J2,—J2 + 1,---,J2 — 1,.72 onto the corresponding directions determined by 
two pairs of Euler angles, (t>i,9i and 4)2, 02- The information contained in the 
tomogram with a dependence on the matrix U of such a form is sufficient to 
reconstruct the density matrix p. But we define the tomogram by Eq.Q to 
use the redundant information on the quantum state of bipartite systems in 
studying the entanglement properties of the system states. We remark that in 
Eq.([3]) we could also use the full unitary group instead of SU{2) and this we 
will do sometimes in the following. 

The tomographic probability distribution of a qudit state ^{m, U) can be 
considered as a column vector yV{U) with components 

WiiU) = W(j, [/), W2(C/) = >V(j - 1, [/), . . . , W2,+i (U) = W(-j, U). (6) 

Since all the components are nonnegative and the normalization condition ([4]) 
holds, from a geometrical point of view the components {Wfe} of the tomo- 
graphic probability vector determine the coordinates {xk} of points belonging 
to a simplex. For a qubit such a simplex is the segment {xi + X2 ~ 1;0 < 



2^1,3^2 < 1} in the plane xi,X2. For a generic qudit the simplex is a polyhedron 
in a {2j + 1)— dimensional space determined by equations: 

23 + 1 

^Xk = l;0<xi,X2,-.-,X2j+i<l. (7) 

fe=i 

Thus, the spin tomogram is a function of a unitary group element U with values 
in a simplex. The linear maps of probability vectors 

W'{U) = MW{U) (8) 

are expressed in terms of matrices M which are known to form the semigroup 
of stochastic matrices. A stochastic matrix is a matrix with nonncgativc entries 
such that the sum of the elements of each column is one. If in addition the sum of 
the elements of each row is one, the matrix is bistochastic. Also the bistochastic 
matrices form a semigroup. In the next section we study the properties of the 
tomograms and stochastic and bistochastic maps on the example of qutrit states. 

3 Qutrit density matrix in tomographic proba- 
bility representation 

To present the semigroup approach to describe generic qudit states we start 
from qutrit states (spin-1 states). The density matrix of a qutrit state is a 
nonncgativc hcrmitian 3x3 matrix with unit trace and it can be presented in 
the product form 

P = U^pUl (9) 

where the 3x3 matrix p is diagonal: 

/5i 

P2 I (10) 

P3 

with nonnegative eigenvalues pk^k ~ 1,2,3, satisfying the normalization condi- 
tion 

3 

E-Sfc^i- (11) 

fe=i 

The columns of the 3x3 unitary matrix Uq are the components of the normalized 
eigenvectors of the density matrix p, i.e. 

Uii Ui2 Ul3 
C/o = I U21 "22 "23 I =||W1,M2,W3| 
"31 "32 "33 



such that 



pUl = piUl , pU2 = P2U2 , PUZ = P3M3 . (12) 



The eigenvalues pk can be considered as the components of a probability vector 



Pi 

P2 
P3 



(13) 



corresponding to a point on the triangle with vertices in (1, 0, 0), (0, 1, 0), (0, 0, 1), 
which is the two dimensional simplex given by Eg. pip with the constraints 
0<Pfc<l, k = 1,2,3. The triangle (qutrit simplex) is embedded in the plane 
determined by Eg. pip , where the coordinates pk arc arbitrary real numbers. 
The linear maps of the qutrit probability vector p are determined by the 3x3 
stochastic matrices M, i.e. 

p*-^^'-Af| (14) 

where the matrix M has real nonnegative entries and the sum of the elements 
of each column is unity. The product of two stochastic matrices is again a 
stochastic matrix. Thus stochastic matrices form a semigroup. Also, any convex 
sum of stochastic matrices is a stochastic matrix. But one can check that the 
inverse matrix AI~^ if it exists is not stochastic because some entries are negative 
numbers. Thus the triangle (qutrit simplex) is invariant under the action of a 
linear stochastic map and the point on the plane pip given by the vector Mp 
belongs to the same triangle. 

The tomogram W{m, U) of the qutrit state with density matrix p given by 
Eq.Q is the probability vector with components 



Wi(C/) 

W2(C/) 
>V3(t/) 



W(+l,C/) = (l UUopU^U^ 
W(0,C/) = /o UUopUlu^ c 
W{-l,U) = l^-l UUqpU^oU'' -1 



(15) 



Direct calculation shows that the above formulae can be written in the following 
form: 

3 

Wk{U) = Y.\{UUo)k^fph (16) 



or in vector form: 

WiU) = Mp 

where the elements of the matrix M read 

A4^ = \{UUo)k„\ 



(17) 



(18) 



As a product of unitary matrices, the matrix UUg is unitary and the sum of 
the (nonnegative) elements of each column and each row of M is one. So, M 
is orthostochastic, a particular bistochastic matrix. The result of Eg.s P?|) . p8p 
means that the qutrit state is determined by a bistochastic map acting on the 
probability 3— vector whose components are the eigenvalues of the density ma- 
trix. The bistochastic matrix M of the map has elements which are the square 



moduli of the elements of a unitary matrix. In turn, the unitary matrix is the 
product of two unitary matrices, one rotating the basis in the space of spin states 
and the other one having columns formed by the eigenvectors of the density ma- 
trix. Prom a geometrical point of view, this means that qutrit states are the 
orbit of the unitary group acting on the points of the triangle (qutrit simplex), 
but the action of the unitary group is made via the action of an orthostochastic 
map. 

In the case of two qubits that construction yields the description of the 
quantum states of the corresponding two spin— 1/2 system by the tomographic 
probability 4— vector VV(C/), where the the unitary matrix U belongs to the 
group C/(4). The spin tomographic probability vector yV{U) has the following 
components 



WiiU) 



W(^,^,C/) = /ii|c/C/o/5C/otf/^ 



1 1 
22 



(19) 



wil,-l,u)Jl^\uUopU^,u^ 






-1 1 



2 2 
-1 -1 
2~^ 



UUopWoU^ 



1 -1 
2~ 
-11 
"2~2 



C/C/opC/gtc/t 



-1-1 



Here the columns of the unitary 4 x 4— matrix Uq are the eigenvectors and p is 
the diagonal form of the density matrix p of the two qubits state. One can easily 
check that the tomogram >V(C/) can be written in the form of the Eq.s P7|) . (|18p 
if the eigenvalues of p are organized as a 4— vector p. This vector belongs to 
a 3— dimensional simplex, a polyhedron, in a 4— dimensional space. The points 
in the simplex belong to the orbit of the group C/(4) in its Cartan subalgebra 
whose elements arc labelled by nonnegative numbers. 



4 The relation to general linear and inhomoge- 
neous general linear groups 

In this section we find a relation of stochastic and bistochastic maps of the tomo- 
graphic probability vectors to Lie groups. If we restrict to invertible stochastic 
and bistochastic matrices and leave out the nonnegativity of their entries, we ob- 
tain a group since the invertible matrices have inverse matrices of the same kind. 
To prove this let us formulate the stochasticity property in terms of stability of 
a vector eb with all its n components equal to 1. The (column) stochasticity 
property of an n x n— matrix M results by the requirement: 



M 



(20) 



where ej = (1, 1, . . . , 1) is the transpose of cq- This requirement is equivalent 
to demand that the sum of the elements of any column of A/ is 1. The identity 



matrix / obviously satisfies the above equation. Since / = MM ^ one has 

ej = ejAf A/-1 = e^M''^ (21) 

and this proves that Eq. (|^D|) is satisfied also by the inverse matrix Af~^, when 
it exists. 

The stochastic matrix M is bistochastic if in addition to Eq. ([^0]) it satisfies 
the condition 

Meo = eo . (22) 

The same previous argument shows that also the inverse matrix M^^ satisfies 
Ea. (|^^ . when it exists. 

The invertible matrices M satisfying equations (PO]) and ([2^ with real en- 
tries of arbitrary sign form a Lie group Gbs- The subset of Qbs of invertible 
matrices M with real nonnegative entries is an open dense subsemigroup of the 
bistochastic matrices semigroup. 

Analogously, the invertible matrices M satisfying equation (PO)) with real 
entries of arbitrary sign form a Lie group Qs -The subset of Gs of invertible 
stochastic matrices is an open dense subsemigroup of the stochastic matrices 
semigroup. 

It is obvious that one can also consider complex matrices, rather than real 
ones, satisfying the above equations. Then one obtains other Lie groups. 

It is clear from Eg. ip^ that the Lie group Gbs is GL{n — 1,M) for the 
n— dimensional case, and GL{2,M.) for the qutrit case. 

Let us discuss the qutrit case in detail to recognize Gbs as Gi(2,M). We 
choose a rotation O, OO = O'^O = / , in the 3— dimensional space such that 



> 

vf 



" 1 ■ 




" ' 


1 


= 





1 




1 



(23) 



For instance, we choose: 



O 



1 


-1 


/2 


V2 


1 


1 


? 


? 



V \/3 \/3 



(24) 



Let us take the bistochastic matrix M of the form 

xi yi zi 

M = \ X2 2/2 Z2 

xi, 2/3 Z'i 



with 
and 



E.^'=-E,j''=-E,^'= = i' 



Xk + Vk + Zk = I \ fc = 1,2,3 



(25) 

(26) 
(27) 



One can check that 

/ a b \ 

M = OMO'^ = c d (28) 

\o l] 
where 

a = -^{xi - X2 - yi + y2) ; (29) 

= (Z2 - zi) ; 



2V2 

^ = 2*-^^^ ~ 2:3-2/3) . 

The 4— parameter group of Eq. ((28|) is the noncompact Lie group GL{2,M.). As 
the determinant of the orthogonal matrix Ois nonzero, the set of matrices M is 
isomorphic to the group of bistochastic matrices M. 

Besides, if one considers the stochastic matrices Lie group Qg, that is ma- 
trices M of the form ([^5]) satisfying Eq.f^H) but not Eg. ipT)) . one can check 

that 

A B m 

M = OMO^ ^\ C D n \ (30) 

1 



where 



A = 7^{xi - X2 - yi + y2) ] (31) 

B = — =(a;i +yi - 2zi -2:2 -y2 + 2Z2) ; 

C = ^(j/3--3); 

D = ^(2z3 -X3 -ya) ; 

1 / 

m = —j=(xi+yi+zi-X2-y2-Z2); 

1 , 
n = ^(1-0:3-^3-23)) ■ 



or equivalently 



Xi 

2/1 



= ^ (3-4 + V2B + V2C + D + V6m + V2n + 2) , (32) 

= - (-3A - V2B + V2C + D- V6m + V2n + 2 

= J (-3A + V2B - V2C + D + \/6m + V2?i + 2 

2/2 = - (SA - V2B - V2C + D- V6m + V2n 
6 V 

zi = - (-2V2B -2D + \/6m + V2n + 
6 V 

Z2 = - (2V2B -2D- V6m + V2n + 2 

From previous formulae one can readily see that pure translations { A = D = 1 
and B — C — ) and pure homogeneous linear transformations [ni = n = 0) 
do not correspond to any stochastic matrix. 

The group of matrices of Eg. ((30)) is just IGL{2, E), isomorphic to the direct 
product of dilation group times ISL{2,R) , which in this particular dimension 
is isomorphic with the Poincare group in two plus one dimensions. In this way 
we have established that the stochastic invertible semigroup of 3 x 3— matrices 
is in one-to-one correspondence with a subset of the group IGL{2,M.). 

Finally, it is clear that this argument holds true for arbitrary qudits (arbi- 
trary spin— j states). 

Now we are able to describe briefly the action of the stochastic and bis- 
tochastic groups on the carrier space M" of their above representation. The 
fiducial vector ej determines an invariant foliation given by 

e^{Mv) = ej(t7) = const ; v e R" (33) 

so that the simplex of probability vectors belong to the leaf ej(w) = 1. Moreover 
that simplex is invariant under the action of the whole semigroup of stochastic 
matrices, and we recall that the orbit of the semigroup of stochastic or even 
bistochastic matrices starting from a vertex is the whole simplex. While the 
action of the stochastic semigroup on the simplex is transitive, the orbit of the 
bisthocastic semigroup starting from a point v — {vi, ...,Vn) is the convex hull 
of all the vectors v-^ whose components are some permutation (w^^, ...,v^^) of 
the components of the given vector v. 

In particular, any point of the simplex can be also reached from a given 
vertex with an invertible stochastic matrix. The same holds true for invertible 
bistochastic matrix, if some lower dimensional set of points of the simplex is 
dropped, in any case the point eo/n cannot be reached. 

From the point of view of matrix analysis (for a general reference of the 
following discussion see, e.g., [H]), a generic stochastic matrix M is a positive 
matrix, i.e. all its entries are positive. Then, by Perron's theorem, we know that 
its spectral radius, which is 1, is an algebraic simple eigenvalue of maximum 



10 



modulus (Perron root of the matrix), with an eigenvector p (Perron vector) 
which may be chosen to be positive, i.e. a probability vector. 

Afp = p ; p, >OVi = l,..,n;y p, = \ (34) 

The fiducial vector eb is just the left Perron vector of the stochastic matrix M 
belonging to the same eigenvalue 1, and the following limit docs exist: 

lim M'= = L := pej ^ \\p,p, . . . ,p\\ . (35) 

k — ^CXD 

In other words, the matrix M^ approaches a limit which is a rank-one stochastic 
matrix whose columns arc the Perron vector of M. As a result, because Lv = p 
for any vector of the leaf e^{v) = 1, we have that 

lim M^v^p (36) 

k — *oo 

independently of the starting point v. In particular, if M is bistochastic, p = 
eo/n. 

We do not insist on the possibility of characterizing that limit in terms of a 
stochastic Markov process. 

When the stochastic matrix M is nonnegative, rather than positive, again a 
nonnegative probability eigenvector p of the eigenvalue 1 exists, and we recover 
the same limit L if the spectral radius of M is the only eigenvalue of maximum 
modulus. 

More generally, if the matrix I — M + L is invertible, we get the same limit 
L by a Cesaro summation: 



N 



1 ^ 

hm -yA/^ = L:=peJ = ||p,p,...,pl|. (37) 



k=l 



That hypothesis is crucial as the following example shows: consider the 
permutation bistochastic matrix 




M= I 1 I (38) 

it gives a Cesaro limit 

. N ( h h ^ 

^--^fc Vo 1 

which does not have the form of L. In any case the columns of the Cesaro limit, 
when it exists, are eigenvectors of M belonging to the eigenvalue 1. 



11 



5 Positive maps 

Having discussed the transformations of probability vectors we are now in the 
position of describing transformations of density states by means of the tomo- 
graphic representation of them. 

The problem of constructing positive maps has a long history and the dy- 
namical maps providing a density matrix which could appear in a process of 
quantum evolution was studied in [5] , see also [T^l US] • 

For finite level systems the space of states is a convex body in the dual space 
of the infinitesimal generators of the unitary group, say u*(7i). For closed quan- 
tum systems the dynamics is conventionally represented by transformations as- 
sociated with one-parameter groups of unitary transformations in Hilbert spaces, 
this formalism however is not appropriate to deal with irreversible behaviours 
like decay of metastable particles, approach to thermodynamic equilibrium and 
more generally to deal with decoherence phaenomena. 

For finite level quantum systems, probability measures are fully described by 
probability vectors, say n— vectors v with nonnegative components whose sum 
is 1: 

Iv: w,„ > Vto = 1, . . . , n; ^^^ Vm = l| ■ (39) 

For those systems, states are identified with n x n— density matrices p which are 
positive semi-definite and satisfy 

p^ = p, Trp = 1 (40) 

Dynamical maps would assign to each p another density state p{t) = M(i)po, 
with M some n^ x n^— matrix satisfying appropriate constraints to guarantee 
that p{t) is still a density state. A dynamical map on states induces a dynamical 
map on probability vectors which may be described by stochastic matrices. 
These maps, on states are called quantum stochastic maps. 

Our aim is to relate stochastic maps acting on probability vectors to quantum 
stochastic maps acting on density states, by using the tomographic setting. This 
approach is different from the standard one |14l[T51ll6| which describes quantum 
stochastic maps as projection of isometries. 

We may parametrize density states by a pair 

P^iUoJ) (41) 

where Uq, as a unitary matrix, provides us with the eigenvectors of p while the 
components of p are the corresponding eigenvalues. Symbolically we could write 

pUo = pUo (42) 

meaning that p acting on the fc— th column of Uq will provide the same col- 
umn multiplied by the corresponding eigenvalue pk- Equivalently, the columns 
of Uq are an orthonormal frame in the Hilbert space and therefore define a fam- 
ily of orthonormal projectors while p gives the weights to attribute to each 
corresponding rank-one projector of the decomposition of p. 
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It should be clear from this description that Uq is determined up to ([/(l))"", 

P=^ Pm \e„i) (em\ (43) 

771—1 

in a given chosen basis of orthonornial vectors. However if the eigenvalues are 
degenerate, the ambiguity increases, replacing each U{1) subgroup with a U{k) 
subgroup depending on the degeneracy of each eigenvalue. Thus, in order to 
define unambiguously the density matrix map of Eq. (j4ip we have to choose a 
"gauge" by fixing the phase factors of Uq and an ordering of both the components 
of p and the columns of Uq so that Eg. ([42]) holds true. 

This particular parametrization allows to deal more easily with quantum 
stochastic maps p ^ p' parametrized in terms of unitary maps and stochas- 
tic maps, providing therefore a different parametrization with respect to those 
normally used in the literature. The density matrix p is mapped onto another 
density matrix p' in such a way that the map is convex linear. This means that 
if two density states are mapped onto other density states 

Pi ^ p'l ; P2 ^ P2 , (44) 

any their convex sum is mapped onto the same convex sum of the images 

Aipi + A2P2 -^ Aip'i + A2P2 . (45) 

In the tomographic framework one has the tomographic map of the density 
matrix pi onto the probability vector of the qudit state 

Pi^>Vi([/) = |t/C/oi|'pi (46) 

and of the density density matrix p2 onto another probability vector of the qudit 
state 

P2 ^ m{U) = \UUo2\^ P2 (47) 

We introduce positive maps of density matrices, parametrized by an unitary 
matrix V and a stochastic matrix M, using this tomographic setting in the 
following way: 

|C/f/oi|'pi ^ |C/C/^i|'p; = |C/y[/oi|'M^i , (48) 

IC/C/02I' P2 ^ |f/C/o2l' P2 = \UVUo2\'' MP2 . 

Thus our map is equivalent to the map of the pair ([/g, p) onto another analogous 
pair (C/q,/5') obtained by acting with an unitary matrix V and a stochastic 
matrix M as: 

{UoJ)^{Ui p') - {VU^ , MP) . (49) 

Since any density matrix p is completely determined by the corresponding pair 
{Uo,p) and as the left actions on the unitary group and on the simplex by 
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stochastic maps are transitive, we have described a class of positive maps which 
act transitively on density matrices. 

To illustrate that the above map is convex linear we discuss in detail the 
simple qubit case. To do this first we find explicit formulae for the maps of 
density matrix p and pair (Uq, p) 



P^iUoJ) 



(50) 



and viceversa 

{UoJ)^p. (51) 

As we have already remarked the map of Eq. ((50| requires a "gauge" choice 
because we are going from a three dimensional manifold to a five dimensional 
one. 

Let the qubit density matrix be 



Then its eigenvalues read 



Pi 



^" P'A , 0<dctp<i. 
P21 P22J 4 



- + -v/l-4detp 



1 1 . 

P2 = ---Vl-4detp. 



2 2 
The eigenvector corresponding to pi is 

\U21 



Ul = 



where 



Mil = 



W21 



2p 



12 



1 + VI - 4dctp - 2/911 
2pi2 



Vi 



yi 



1 + vT^^ldetp - 2pii 
The eigenvector corresponding to p2 is 

A^12 



"2 = 



\U22 



where 



W12 



U22 



2p 



12 



1 - VI - 4dctp - 2/911 

2pi2 



2/2 



y2 



1 - VI - 4dctp - 2/911 



(52) 

(53) 
(54) 

(55) 

(56) 
(57) 

(58) 

(59) 
(60) 
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and the phases have been chosen to be zero. Thus given the qubit matrix p of 
Eq. (j52p one has the corresponding pair of unitary matrix 



^o = r" "^M (61) 

\U21 U22J ^ ' 

([55)) . ([50)1 and the point on simplex given by the 

'.- (J) (e^) 

One can check that the matrix p of Eq. ([5^ has 



with components ([5 
probabihty vector 



with components ([2 
the representation 



P = f/o 



Pi 
h 



ul 



(63) 



Let us consider the inverse problem. Namely, given a unitary matrix Vb of 
the form (pTj) and a probability vector 



X 



, Xl + X2 = 1 



describing a point on the simplex, let us construct the matrix 



x = K, 



Xl 







X2 

whose matrix elements explicitly are: 



V^ 



Xii X12 
X21 X22 



Xll = l^lll Xl + I^12| X2 = l-X22; 

X12 = ^"11^2*1X1 + ^"12^/72X2 = X21 ■ 



(64) 



(65) 



(66) 

(67) 



Thus, given the pair (Vo,x), one has a density matrix x with the previous 
matrix elements. One can see that the map p — *■ {Uo,p) is nonlinear. The 
formulae ([55]) , (p7)) provide a polynomial dependence of the elements of density 
matrix on unitary matrix and probability vector components. The dependence 
of density matrix elements on probability vector components is linear and on 
unitary matrix elements Vjk is quadratic. It is obvious that the inverse transform 
p -* {Uo,p) is also nonlinear as formulae ^, dM]), dSHl), dSZ]), dSSD, dSO]) show. 
Moreover formula (|66p shows how two U(l) elements disappear in defining Xjk- 

Let us discuss now the representation of a convex sum of two density matrices 
p^^' and p"' given by 



p = Aip(i) + X2P^^\ Ai + A2 = 1, < Aj < 1 
The matrix elements of this matrix read 

Pjk = >^ipfk + >^2pfJ 



(68) 
(69) 
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The two eigenvectors of this matrix are 

pi,2 = i ± iy^l - 4A2dctp(i) - 4A2dctp(2) - 4A1A27 
where 

(1) (2) , (2) (1) (1) (2) (2) (1) 

7 = P11V22 + Pii P22 - P12P21 - P12P21 
One can check the following equality of probability vectors: 



C/C/o(A) p(A) = Ai t/C/^'^ ' pi + A2 UU^^^ ' p2 



(70) 
(71) 

(72) 



where p{\) is the probability vector with components given by Eq.s ((70|) . ([71]) . 
The mritary matrix Uq{\) is obtained by formula ([6T|) with entries ([56]) . ([57]) . 
([59]) . ([60]) after the replacements 



P12 - 

Pii - 
detp 



^lPl2 



A2P12 



(73) 



^iPll + ^2Pii 

-* Ajdetp^^) + A^detp(2) + A1A27 



where 7 is given by Eq. (j7ip . Thus we checked directly the convex superposi- 
tion formula for two tomographic probability vectors corresponding to a convex 
mixture of two density matrix p^^' and p^"^' for two qubits. From Ea. (|72p one 
can readily see that 



p*(A) = Ai U^^XM^^ 



Pi + A2 



C/o-i(A)C/f 



P2- 



(74) 



Thus we have illustrated in detail that the map of pairs {Uo,p) — > {VUo,Mp) 
by means of left action on unitary group and stochastic matrix map on simplex 
is convex linear. It is worth to note that the map (Uo,p) —i- (Uq,p) , which 
is just the transposition of the density matrix map, is not equivalent to left 
multiplication by a matrix. So, we can extend the map using automorphism 
groups of both the unitary matrices and simplex points. 



6 Entanglement 

Stochastic maps may be used successfully to characterize the entanglement of 
states for a composite physical system in the tomographic scheme as the follow- 
ing example shows. 

To study entangled states we consider two qubit states with state vector 



w 



1 

V2 



The density matrix of this state reads 



I) 


+ 


-I) 


1 

+2 


/o 






1 


0' 

1 







1 


1 


7 


Vo 





0, 





(75) 



(76) 
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so its four eigenvalues yield the probability vector 














p = 


1 





7 




while the corresponding eigenvectors form Uq : 




f/o = 


/I 

Vo 


- 



1 

"71 



0\ 





1/ 



(77) 



(78) 



In view of this, the tomographic probability vector of the two qubit state has 
the following form: 



W{U) 



\Ul2 + "131 
\U22 +U23P 
1^32 +W33P 
\ui2 +U43P 



(79) 



where the u's are the matrix elements of U. One can recognize the entanglement 
of the state with this tomogram, calculating the stochastic matrix M whose 
columns are four probability vectors VV(C/ab), VV(C/ac), VV(C/db), VV(C/(ic), where 
each unitary matrix Uhk is a tensor product of two 2x2 unitary matrices 

Uh ® Uk ■■ 

Uhk ^Uh®Uk , (h^a.d; k^b,c) . (80) 



Eventually, the stochastic matrix M has the form 

Xdb 
2 '^db 
2 '^db 

Xdb 



M 



1 Xab 


^ac 


2 ^ab 
2 -^ab 


2 "^ac 
2 -^ac 


\ Xab 


^ac 



2 ~^ •^dc 
2 •^dc 

•^dc 



) 



(81) 



where xuk is a trigonometric function of the Euler angles determining Uh,Uk ■ 
To evaluate the Bell number B satisfying Bell's inequality [6l [7] 



B <2 



one has to calculate [17] a trace: 



where E is the following matrix 



E 



B- 


= Tr(j 


VIE) 




x: 

/I 

1 


-1 
-1 


-1 
-1 


1\ 

1 


1 


-1 


-1 


1 


V-1 


1 


1 


-V 



(82) 



(83) 
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The elements of the matrix M are functions of four directions. The Cirelson 
bound[8] for the BeU number is 2-\/2 and can be achieved in entangled states 
only. One can check that the bound is achieved when 

In that case one has 

2x/2 = 4(4x-l) ; x = ^^^ . (85) 

o 

Thus the universal stochastic matrix corresponding to an arbitrary, maximally 
entangled state of two qubit is 



M 



/ 2+x/2 2+v^ 2+\/2 2-\/2 \ 
/ 8 8 8 8 

2-\/2 2-\/2 2-\/2 2+V2 

8 8 8 8 

2-\/2 2-\/2 2-\/2 2+V2 

8 8 8 8 

\ 2+^2 2+^/2 2+\/2 2-\/2 / 

\ 8 8 8 8 / 



(86) 



Thus we have clarified what is the relation of stochastic matrices with Bell 
inequality violation for entangled states of two qubits. 

7 Conclusions 

Wc point out the main results of this paper. We have established a connection 
of spin tomograms with stochastic maps acting on a simplex and unitary group 
elements. We have shown that stochastic and bistochastic n x n matrices have 
a dense intersection with the Lie groups IGL (n — 1 , R) and GL{n — 1 , M) respec- 
tively. We have constructed positive maps of density states as maps determined 
by pairs of unitary matrices and stochastic matrices. We have demonstrated 
that for entangled two qubit states the Cirelson bound for the Bell number is 
associated with a universal stochastic matrix. 

To conclude, we observe that in description of quantum states an important 
role is played by unitary irreducible representations of Lie groups. One can 
see that starting with unitary representations of the GL{2,M.) group which are 
infinite-dimensional and restricting the representations to the semigroup subset 
one obtains unitary representations of invertiblc bistochastic matrix semigroup. 
Analogously the unitary representations of the stochastic semigroup can be ob- 
tained by restricting to the semigroup subset the infinite-dimensional unitary 
representations of the inhomogeneous linear group (cither real or complex) . The 
problem of irreducibility or rcducibility of the semigroup representations has to 
be analyzed further. 

Finally, we hope to clarify the relations between the constructed positive 
map and other existing approaches in a future work. 
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